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ABSTRACT 

Time-delay autosynchronization (TDAS) can be used to stabilize unstable periodic orbits in 
dynamical systems. The technique involves continuous feedback of signals delayed by the orbit's 
period. One variant, ETDAS, uses information further in the past. In both cases, the feedback 
signal vanishes on the target periodic orbit and hence the stabilized periodic orbit is one of the 
original dynamical system. Furthermore, this control method only requires the knowledge of the 
period of the unstable orbit. 

In general, the amount of feedback gain needed to achieve stabilization varies with the 
bifurcation parameter(s) of the system, resulting in a domain of control which can be computed 
without having to deal with the explicit integration of time-delay equations. 

In this paper we compute the domain of control of the unstable periodic orbits of the buck 
converter. The simplicity of the nonlinearity of this converter allows us to obtain a closed 
analytical expression for the curve g : S 1 — > C whose index determines the stability. We perform 
detailed numerical explorations of this index for both period-1 and period-2 orbits for several 
configurations of the feedback scheme. We run several simulations of the controlled system and 
discuss the results. 



1 Introduction 



Control of chaos, meaning supression of the chaotic regime in a system by means of a small, 
time-dependent perturbation, has been a subject of interest in recent years ^2]- In [0], it was 
pointed out that the many unstable periodic orbits (UPOs) embedded in a strange attractor 
could be used to produce regular behaviour to the advantage of engineers trying to control 
nonlinear systems in which chaotic fluctuations are present but undesirable. It has also been 
speculated that UPOs may play a key role in the regulation of complex biological systems. 

Two main groups of methods of control of chaos such that the feedback perturbation vanishes 
on the target orbit have been considered in the literature: 

• The method proposed by Ott, Grebogi and Yorke where small perturbations to an 
accessible parameter are introduced. The method exploits the fact that during its wander- 
ing over the strange attractor, the system will eventually come near the target UPO on a 
given Poincare section. When this happens, and only then, a small perturbation is applied 
to the parameter so as to make the orbit land on the stable variety of the target orbit the 
next time it crosses the Poincare section. As a drawback, the method is not suitable for 
complex systems, since a nontrivial computer analysis must be performed at each crossing 
of the Poincare section. Also, small noise can drive the orbit away from the target orbit, 
and the control method must then wait for a while until the system comes near to the 
target orbit again. 

• The method proposed by Pyragas ^0] > called time-delayed autosynchronization (TDAS), 
involves a control signal formed with the difference between the current state of the system 
and the state of the system delayed by one period of the UPO. One variant, ETDAS, 
proposed by Socolar et al ^3], uses a particular linear combination of signals from the 
system delayed by integer multiples of the UPO's period. Still another variant, proposed 
by de Sousa Vieira et al ^1] uses a nonlinear function of the difference between the 
present state and the delayed state. TDAS and its variants have the advantage that 
the only information needed about the target orbit is its period and that no computer 
processing must be done to generate the control signal. The method has even been applied 
to systems described by partial differential equations |2j. In general, the feedback gain 
which succesfully stabilizes the orbit lies in a finite, and often narrow, orbit-dependent 
range. In the space of the feedback gain and the bifurcation parameter(s) of the system, 
the region where the TDAS can be applied with success is called the domain of control. 
In a method was proposed to compute the domain of control of a given system without 
having to explicitly integrate the resulting time-delay equations, which is a nontrivial 
matter due to the choice of initial conditions 8 . Essentially, the method reduces to the 
computation of the index around the origin of a curve in the complex plane. 

In this paper we address the problem of stabilizing the UPOs of the buck converter by 
means of ETDAS. The chaotic behaviour of this converter and the computation of its UPOs 
have been extensively studied in [3], [H] and [7j, and methods of control based in OGY techniques 
and others have been analyzed [I][E]- The main result of our work is that the function g fron 
the unit circle to the complex plane whose index determines the success of ETDAS can be 
analytically computed for the buck converter. The index can then be numerically evaluated and 
the domain of control can be easily constructed. 

The paper is organized as follows. In Section |2] we present the three schemes that we have 
studied to feed back the control signal to the buck converter. In Section we review the method 
os Bleich and Socolar to compute the domain of control and specialize it to variable structure 
systems. In Section |I] we analytically compute g for the period-1 orbits of the buck converter 
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(a) First scheme (b) Second scheme (c) Third scheme 

Figure 1: The three feedback schemes for time-delay autosyncronization of the buck converter 

and the three proposed schemes. In Section El we perform numerical explorations of the exact 
g functions obtained in Section |1] and draw the corresponding domains of control. In Section 
1^1 we partially extend our results to higher period orbits of the buck converter and treat in 
some detail the 2-periodic orbits. In Section we perform explicit simulations of the ETDAS 
controlled buck and get numerical, independent confirmation of our analytical results, and at the 
same time illustrate the effectiveness of the ETDAS technique. Finally, we state our conclusions 
in Section ISJ 



2 Time-delayed feedback for the buck converter 

Figure n shows the basic scheme of the PWM controlled buck converter together with the three 
different schemes that we have studied to implement ETDAS. v(t) is compared with a periodic 
ramp 

r (t) = V ief + — + ^ t, 

with t evaluated mod T, and the switch S is opened if v(t) > r(t) and closed if v{t) < r(t). If 
we assume that the inductor current i(t) is always positive, the system has only two topologies 
and is described by 

dv 1 1 di IE., 

— = v H 1, — = v H ult) 

dt RC C ' dt L L 

where u(t) = 1 — 9{v{t) — r(t)). The input voltage E acts as a bifurcation parameter. Notice 
that the system is nonlinear only because of the changing topology. We refer to jSj, [HI and [7] 
for a detailed description of the bifurcation diagramm of this system and the chaotic regime. In 
all the numerical simulations we will use the values R = 22 0,, C = 47 fiF, L = 20 mH, a = 8.4, 
T = 400 fjs, V L = 3.8 V, Vu = 8.2 V, V Tei = 11.3 V and E varying in the range [20, 35] V. 
The control signal Av(t) is given by 



Av(t) = j] v(t) -(l-r)^2 r k ^ X v{t - kr) 



\ k=l / 

where r\ is a (dimensionless) feedback gain and r € [0, 1) determines the relative weight of the 
increasingly delayed contributions. The case r = corresponds to TDAS. r > is the period 
of the target UPO which, in our case, is always a multiple of the ramp period T. Notice that 
Av(t) = if v(t) is periodic with period r. 

We will be primarily concerned with the stabilization of UPOs which cross the ramp exactly 
once every period, and such that at the beginning of every period one has v{t) > r{t). More 
general situations can be easily treated as well, as will become clear in Section 0J 
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3 ETDAS for variable structure systems 



The PWM controlled buck converter described in Section as well as other PWM controlled 
DC-DC converters, can be considered as particular cases of dynamical systems with equations 
of the form 

x(t)=A(t,[x])x(t) + b(t,[x}), (1) 

where x,b € M. n , A <G M(lR n ). We use the notations A(t, [x]) and b(t, [x]) to indicate that both 
A and b are local functionals of x. For the PWM controlled converters, A and b are piecewise 
constant, tipically changing their values when a linear function of x crosses a given periodic 
function of t. 

Consider now a T— periodic orbit x*(t) of this system and a nearby orbit x(t). We want to 
study the evolution of y(t) = x(t) — x*(t). We will be mainly interested in unstable orbits x*(t), 
and our goal will be to modify the right-hand side of (JIJ) so as to render x*(t) stable, that is 
rimt_ > _|_ 0O y(t) = for x{t) initially close enough to x*(i). We will do so by means of extended 
time-delay autosynchronization and consider the following equation: 

+oo 

x(t) = A(t, [x])x{t) + b(t, [x]) + r)M(t, [x])(x(t) - (1 - r) ^ r fc_1 x(t - jfer)), (2) 

k=l 

where M € M(M n ) is a matrix indicating how the delayed signal is feeded back to the system 
and 7] is the strength of the feedback gain. 

We will study the evolution of y(t) to first order in y under @. One has 

y(t) = A(t,[x])x(t)-A(t,[x*})x*(t) + b(t,[x])-b(t,[x*}) 

+oo 

+ V M(t, [x*])(y(t) - (1 - r) £V-yt - fer)) + 0(y 2 ) 

k=l 

where the r— periodicity of x* has been used and we have thrown away terms of higher order in 
y. Expanding the functionals around x* gives 

ij{t) = A(t,[x*])y(t) + 

+oo 

+ r ? M(i,[x*])( 2 /(t)-(l-r)^r fc - 1 y(t-A;r)) + 0(y 2 ). (3) 

k=l 

Since the functionals are local, the functional derivatives will yield delta functions and, forget- 
ting about the higher order terms, we get a variable coefficient, linear time-delayed differential 
equation for y which, with suitable definitions, can be written in the form 

+oo 

y(t) = A (t)y(t) + Mt)v(t ~ kr), (4) 

k=l 

where the time dependence of the coefficients is known since x*(t) is known. Notice also that 
these coefficients are periodic functions of time, since their explicit dependence of time is periodic 
and x* is also periodic. From (j2J) and Q) one can easily read 

A k (t) = -(l-r)7 1 M(t,[x*])r k - 1 , k = 1,2, (5) 

We are not interested in the general solution of the time-delayed equation (j3J) , but rather we 
would like to know if its zero solution is assimptotically stable. To this end, we look for solutions 
of the form 

y{t)= Px {t)e xt l\ 



5A{t, [x]) 
6x k (t') 



yk(t') dt' 



x*(t)+ [ 
Jo 



T Sb(t,[x}) 
5x k (t>) 



Vk{t') dt' 
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with A E C and p\(t + r) = p\(t). This yields for p\ an ordinary differential equation 



p x (t) = (A (t) - -I) Px (t) + V e- kX A k (t) Px (t), 

k=l 

whose solution for a given initial condition can be expressed in terms of an evolution operator 
U x (t) defined by 

px(t) = e- xt / T U x (t)p x (0) 

and satisfying the equation 

/ +°° \ 



U x (t) = \ Ao(t) + £ e~ kX A k (t)j U x (t) (6) 

with initial condition U\(Q) = I. The general solution to @ can be formally expressed as 

/ r-t +°° \ 

Ux{t) = Texp f ^ (Ao(t') + J>- fe \4 fe (0) di'J , (7) 

where T stands for time ordered product (this formal solution is also known as Peano-Baker 
series in the mathematical literature jllj : it boils down to the standard exponential matrix if the 
coefficients of the differential equation are constant). In any case, U\{t) retains the fundamental 
properties 

Ux(t 1 +t 2 ) = Ux(t 1 )Ux(t 2 ), U x (-t) = U^(t). (8) 
Using (jHJ) the condition p\(t + r) = p\(t) is easily seen to be equivalent to 

(e- A C/ A (r)-l) p x (0) = 

which implies 

det (e~ x Ux{T) -l\ = (9) 

This equation determines the values of A such that p\(t) is a solution of our equation. As we 
want y(t) = e xt / T px{t) go asymptotically to zero, we must demand that 3?A < for all the 
solutions of ©. Defining the Floquet multiplier fi = e x and t/(r;/x -1 ) = U\(t), equation (jHJ) 
becomes 

^- 1 ) = detOi _1 l7(T;/i- 1 )-I) = (10) 

with 

Ufafi- 1 ) =Texp U (A (t) + Yl^' k Mt)) dtj . (11) 

Summing up, stability of x* is equivalent to the requirement that all the zeros of g(n~ l ) lie 
outside the unit circle (||/^ _1 || > 1 O 3?A < 0). Using © the series in (fTT|) can be summed to 
yield 

U(t; fi- 1 ) = Texp (j\Mt) " (1 " r) V M(t, [x*]) ± * ) dtj . (12) 

For r < 1, U(t; /^ _1 ) as a function of [jl , and hence g , (// _1 ), has no poles inside the unit circle. 
By a well-known theorem of complex analysis, the number of zeros of g(fi~ 1 ) inside the unit 
circle equals the index with respect to the origin of the curve traced by g([i~ 1 ) when fi^ 1 runs 
over the unit circle. Therefore, the solution will be stable if and only if this index is zero. This 
way of computing the number of zeros inside the unit circle is numerically preferred in front of 
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the obvious option of actually computing the zeros, which in fact are infinite (this is due to the 
time-delay character of the original equation). 

The equations derived so far are completely general. In practice, however, analytical com- 
putation of U{r;fi~ l ) using 1)12(1 is not possible except for very special cases, so one must fall 
back to the defining differential equation © and numerically integrate it between and r. For 
the buck converter, the functional derivatives appearing in (J2J) are specially simple (in fact, A is 
constant), and the differential equation © can be analytically integrated. We will do so in the 
next section, for several choices of M. 

4 Analytical computation of g(fi^ 1 ) for the buck converter 

For the buck converter, one has 

A(t, [*]) = ( ~JF f ) = A and b(t, [x]) = ( |) (1 - 9(v(t) - r(t))) 

so that all the functional derivatives of A(t, [x]) are zero and the only nonzero functional deriva- 
tive of b(t, [x]) is 

5b2(t, \x\) £/.... , , 
Then, equation © becomes, with r = T, 

ij{t) = Ay(t) - | ( J ° ) y(t)5(v*(t) - r(t)) + V M(t, [x*]) (y(t) - (1 - r) £ r^V* - fcT)J 

(13) 

It is easy to see that the three feedback schemes described in Section |2] correspond to the 
following choices for the matrix M: 



1. 



2. 



3. 



Mfe [*]) = ( ^ o)^ M i 
M(i, [*]) = ( j| J ) = M 2 



M(i, [*]) = (1 - - r(t)) ( ? ° J = (1 - 0(w(t) - r(i)) M 3 

Integration of the resulting differential equation is very similar for cases 1 and 2, while 3 is 
slightly more involved. We will present the explicit computation of g{^ 1 ) for cases 1 and 3 and 
deduce the result for case 2 from the one obtained in case 1. 
In the first feedback scheme one has 

Mt) = A + ^-j(i J )*(«*(*) -r(t)) 
A k (t) = -^(l-ry^Mx, A; = 1,2,... 
and the differential equation for U(t) is, adding up the geometrical series, 



L L 



S(v*(t)-r(t)) 
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If we write 

ui(t) u 2 {t) 

u 3 (t) U A {t) 

we get two uncoupled systems of dimension two 



U{t) 



«l(t) = 


-a\ui(t) + 






u 2 (t) = 


-a\U2(t) + 






«s(t) = 


— «l(t) - 


j5(v*(t) 


-r(t)) Ul (t) 


u 4 (t) = 


-i«2(t) " 


-*(«*(<) 


-r(t))u 2 (t) 



where 



a i 



1 

RC 



— — 1 — r/ + rjfi 



,-i 



with the initial conditions ui(0) 
to solve twice the single system 

x(i) 

m 



1 — r/i 

1, 1*2(0) = 0, 2/3(0) = 0, 2x4(0) = 1. Therefore, we only need 



1 

-a x x(t) + — y{t) 
~x{t)--8{v*(t)-r{t))x{t) 



(15) 



for t € [0, T] with the two sets of initial conditions (1, 0) and (0, 1). If we assume that for each 
cycle of the auxiliary signal r(t) there is one and only one t = t c E (0, T) such that the system 
switches topology, then the delta function appearing in the above equation is 



5(v*(t)-r(t)) = (36(t-t c ), 



where 
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\v*(t c )-r{t c )\ 

is the inverse of the absolute value of the slope with which v*(t) crosses the ramp. This is well 
defined because v(t), in contrast to i(t), which is not differentiable at t = t c , is everywhere 
differentiable. We get thus the time-varying linear system 



±{t) 

m 



-aix{t) + ^y(t) 



L 



x(t) 



0E_ 



5(t - t c )x(t) 



(16) 
(17) 



Although this is a variable coefficient system, it can be solved by Laplace transform since the 
delta function makes it trivial to compute the transform of the last term of (J17|). We skip the 
details and write down the general solution for t £ [0, T]: 



x(t) 



y(t) 



-x(0)^-^ + y(0) 
271 71 



e -|K+7i)t 



E(3 



■x{t c )6{t - t c 



e -!(a l -7l)(*--t c ) _ e -|(ai+7i)(t-tc) 



-*(o)-^ + y( o) ai + 71 



L71 



271 



e 2 



|(ai-7i)t 



+ x(0) 



L71 



1/(0) 



«i ~ 7i 
271 



e 2 



|(ai+7i)t 



E(3 
2Z7" 



-s(tc)0(t - i c ) ((ai + 7l ) e -|^-^)(*-*=) _ ( ai _ 7l ) e -| 
1 v 



|(ai+ 7 i)(t-tc) 
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where 



7i 



LC 



and x(t c ) is denned evaluating the expression for x(t) at t = t c . Notice that x(t) is continuous 
at t = t c , while y(t) has a jump of value —(3E/Lx(t c ), as follows from (|T7|) . 

ow putting x(0) = 1, y(0) = we obtain, at t = T, u\(T) and us(T), while x(0) = and 
y(0) = 1 yield us U2(T) and u^{T). Then we can compute 



det 



//"^(T) M _1 «4(T) - 1 
( W1 (T)«4(T) - U2(r)ti 3 (r)) - /i- 1 («i(T) + n 4 (T)) + 1 



(18) 



which, after a little algebra, yields 



The term that multiplies pT 2 could have been easily computed without solving the system, since 
it equals the Wronskian, W, and satisfies W = — a{W , with W(0) = 1. 

Notice that, although 71 introduces a line cut where its argument vanishes, if we analytically 
extend the above expression by its series expansion, only even powers of 71 do appear and thus 
<?i(ti _1 ) only has pole-type singularities (outside the unit circle for r < 1). 

In case 2, the equation obeyed by U(t) is 



U(t) 



1 1— r 



-ryr 



1 

C 



.i ( l + r? _ w -i_J- )_|^;*(i)-r(t)) o 



tr(t) 



Comparing with case 1, the differences amount to changing the sign of 77 and then replacing 1/L 

1 r 



by 1/^(1 +r] — rjfi 1 y~r^T") without changing the combination £?/L. Thus, the index function 



in this case is 



52 (M *) = M ^ e 



-2-a 2 T 



cosh^T) - ^ sinh(i 72 T) j + 1 (20) 



with 



«2 



RC 



1 + 77 - 77/j 1 



1 — rfj, 1 



and 



72 



4 / ! 1 

•— 1 + 77 - 77/j 1 - 



f 2 LC \" ' '' ' ,r l-rfi-\ 
For the third scheme the evolution operator U{t) obeys the equation 

/ 



U{t) 



1 

RC 



1 

C 



V -T ~ f *(«*(*) " K*)) + £(1 " *(**(«) " r(t)) (l - /^i^r) 



U(t) 



Following the same steps as in the previous cases, we arrive at the time- varying linear system 

±{t) = ~Jic x{t) + c~ y(t) 



\x{t) - ^5{t - t c ) x {t) + ^e(t - t c )(i - 1 -L- ) x (t), 



L 



L 



1 — r/i 
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which we must solve for the initial conditions (1, 0) and (0, 1). For t c < t < T the solution is 

+ ^(t c )^-y(t c )-^)e-|(^)(^) 
y(t) = (-^(t c )-^ + 2/ (t c )^)e-^^3)(^) 



where 

z(i c ) = ^_ x (0)^ + y (0)^ e -4(^+^(0)^-y(0)^)e-i^>*» 

vftO = (-*(o)^ + y(o)^ 2 ) e4(a " 7)tc + (*(o)^ - y(°)^f) e " |(a+7)tc " x x(tc) 

(here y(t c ) really means y(tj)) and 



After some calculations one can get 

E/3 



g 3 (^) = ^e- T - 2^e-^ T (A - ^B) + 1, (21) 



where 



1.x wl , m . xx . W 2 o^ + M . wl . x . , ,1 



i4 = cosh(- 7 t c ) cosh(- 73 (r - i c )) + [a*- 2——^ sinh(- 7 i c ) sinh(- 73 (T - t c )) 

2 2 773 V CLL 3 J 2 2 

B = -sinh(-7i c ) cosh(^7 3 (T - t c )) + — cosh(-7t c ) sinh(-7 3 (T - t c )) 
7 2 2 73 2 2 

The main new feature of this expression with respect to g\ and 52 is the explicit appearance 
of t c . As a check, one can see that <? 3 (/U _1 ) reduces to g\ (^ _1 ) if one sets di = a, 73 = 71 = 7, 
as it must be. 



5 Numerical analysis of g(fi : ) 

We have numerically evaluated the index functions gr(// _1 ) of the previous Section. Figure [21 
shows the domains of control for the two first feedback schemes and r = 0.0, r = 0.6 and 
r = 0.9, and Figure 01 corresponds to the third feedback scheme and the same values of r. The 
black regions are those yielding index 0, when the time-delayed feedback succesfully stabilizes 
the T-periodic orbit. Several features must be highlighted: 

1. There are no big differences between the first and second schemes, apart from the sign of 
the feedback. The domain of control of the first scheme is slightly broader. 

2. Use of extended time-delay feedback does not improve the domain of control for the two 
first schemes. The index 1 zone expands with r at the expense of both the index and 
index 2 zones so, in fact, the stable region diminishes. 
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3. The third scheme is clearly inferior, as the stabilizing region is bounded and for low 
values of r does not reach into the chaotic region. This is not surprising, since in this 
scheme the feedback signal only works during one of the topologies. In this scheme, use 
of extended time-delay feedback broadens the stabilizing region and allows control in the 
chaotic regime. 




(a) First scheme, r — 0.0 (b) First scheme, r = 0.6 (c) First scheme, r = 0.9 




(d) Second scheme, r = 0.0 (e) Second scheme, r — 0.6 (f) Second scheme, r = 0.9 



Figure 2: Domains of control for the two first feedback schemes and several values of r. Vertical 
axis: ij from —10 to 10 . Horizontal axis: E from 20 V to 35 V. Black = (stable), Grey = 1 
(unstable), White = 2 (unstable). 



6 Stabilization of higher order orbits 

In this Section we will show how to generalize the analytical results of Section |I] to more general 
situations than period-1 UPOs with a single crossing per period. We will only work wirh the 
first feedback scheme. 

Consider for instance a period-n orbit with a single crossing per period at times t\ £ (0, T), 
(^2 £ (0, 2T),. . ., t n E ((n — 1)T, nT). Now we have to integrate U(t) over [0, nT] and the system 
changes its topology at T, 2T, . . . , (n — 1)T in addition to t\,t2, ■ ■ ■ ,t n . However, both x*(t) and 
x*(t) + y(t) change the topology of the system at exactly T, 2T, . . . , (n — 1)T, so we may forget 
about these changes since they do not take the orbits apart. Mathematically, using the notation 
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(a) r = 0.0 



(b) r = 0.6 



(c) r = 0.9 



Figure 3: Domains of control for the third feedback scheme and several values of r. Vertical 
axis: r\ from —10 to 50. Horizontal axis: E from 20 V to 35 V. Black = (stable), Grey = 1 
(unstable), White = 2 (unstable). 



of Section|3J this is reflected by the fact that only 5b(t, [x])/5v(t') enters the computations and 
that this functional derivative is zero for all the variations of b(t, [x]) which do not depend on 
v, such as those occurrying at the end of period. 

In the first feedback scheme, we may now proceed to equations (|15f). which we will have to 
integrate between and nT for the two sets of initial conditions (1,0) and (0, 1). For each cycle 
of the auxiliar ramp r{t) we have a zero of the delta function and hence 

n 

5(v*(t)-r(t)) = ^f3 k 5(t-t k ), 

k=l 

where 

Pk= \v*(t k )-r(t k )l 

Integration is now easily performed by means of a Laplace transform and we get the result 

x(t) = f-x(O)^-^ + y(0)-V) e -*<*-*» + Uo)2±2. _ y(0 ) J-) e -¥^)t 
V 2 7i 7iW V 2 7i 7iW 

P n 

' » '■ — '■ ■ * I - k '—n)(t-t k ) _ e -I( 01 +7i)(t-t fe ) 



V(t) = f-x(0)-L + 2/(0)^^ e-K--^)* + ^(o) J_ _ ,(0)^-^) e-l^+TOt 
V 7i 271 y v 71 L 271 y 



£ n 



- (Pkx{t k )9{t - t k ) ((ai + 7i)e"^ (ai 



fc =i 



-7i)(*-*fc) _ ( ai _ 7l ) e -|(ai+7l)(*-*fe) 



with the x(tfc) computed recursively from the expression for rr(i). 

We will present explicit expressions for g(^ _1 ) only for the case n = 2. One gets 

- 2 / u~ 1 e- aiT f cosh( 7 iT) - -^(/?i + ft) sinh( 7l r) 

+ ^^(cosh( 7 iT) - cosh( 7l (t 1 - t 2 + T)))) + 1. (22) 
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Domains of control for this case with r = 0.0, r = 0.6 and r = 0.9 are presented in Figure 0J 
and the same comments done for the period-1 orbits do apply here. Other situations, such as 
period-1 orbits with several crossings per period, can be treated along the same lines developed 
in this Section. 




(a) r = 0.0 (b) r = 0.6 (c) r = 0.9 



Figure 4: Domains of control for the first feedback scheme and period-2 orbits and several values 
of r. Vertical axis: rj from —10 to 10. Horizontal axis: E from 25 V to 35 V. Black = (stable), 
Dark gray = 1 (unstable), Light gray = 2 (unstable), White = 3 (unstable). 



7 ETDAS simulations 

In this Section we report the results of several simulations of the time-delay controlled system. 

Figure 03 shows a typical numerical simulation of the time-delay feedback control method. 
A chaotic orbit of the system (rj = 0) is also shown for reference. The feedback starts to act 
after the first period and stabilizes the orbit in less than 10 periods, i.e., 4 ms for the system 
considered. 

Figure HJa) shows a simulation of an ETDAS with r = 0.6, rj = —5.0 and E = 33V using 
the first scheme. The series in (j2J) has been truncated to the first 80 terms. The figure shows 
the first 100 periods, and the feedback starts to act after t = 80 • T. Notice that the collapse of 
the chaotic regime to the stabilized period-1 orbit is nearly instantaneous. Another ETDAS, this 
time using the third scheme and r = 0.6, r\ = 6.0 and E = 26V is depicted in Figure Efb). The 
system first stabilizes on the period-2 stable orbit. The time-delay control starts to act again 
after t = 0.032s and, after a short transient, the period-1 UPO is stabilized. 

In order to check the frontier between the zones of the domain of control, we have numerically 
integrated the time-delay feedback equations for several values of the bifurcation parameter E 
and the feedback gain rj on both sides of the analitically computed frontier, although some times 
numerical integration errors may produce a wrong result if rj is too close to the frontier (this 
is particularly true of the third feedback scheme, where numerical errors during the v(t) < r(t) 
topology are amplified during the phase v(t) > r(t), when the control does not act). 

One of those checks is presented in Figure 0a), which corresponds to input voltage E = 30 
with the first feedback scheme and r = 0. The solid line corresponds to rj = — 1.3 and the dashed 
one to rj = —1.2. The expression (|19j) for g^ji^ 1 ) predicts index for the former and index 1 
for the later. The time span in the figure corresponds to four periods of the auxiliary ramp and 
the vertical axis represents the capacitor voltage. We see indeed that rj = —1.3 stabilizes the 
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-| o I i i i i i i i i i 

0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 

Figure 5: Load voltage chaotic waveform for E = 35V (solid line) and time-delay feedback 
controlled orbit in the second scheme with r = and rj = 4.0 (dashed line). 50 periods of the 
auxiliary ramp are shown. 





11.8 ' 1 1 1 1 1 ' 1 11.7 1 1 1 1 1 1 1 

0.035 0.0352 0.0354 0.0356 0.0358 0.036 0.0362 0.0364 0.0366 0.06 0.0605 0.061 0.0615 0.062 0.0625 0.063 

(a) (b) 
Figure 7: Simulation checks of the analytically computed limits of the control domains 



system to the period-1 orbit, while 77 = —1.2 does not. In fact, 77 = —1.2 produces a period-2 
orbit which, however, is not the stable period-2 orbit of the uncontrolled system which exists for 
this value of E and is also represented in the figure: there is some kind of competition between 
the stable period-2 orbit and the time-delay feedback (which is not zero on this period-2 orbit) 
which finally neither stabilizes to the unstable period-1 orbit nor falls on the system's stable 
period-2 orbit. 

Another check, this time for period-2 orbits, is illustrated by Figure Efb), corresponding to 
E = 32.5 and r = 0.0. Eight periods of the auxiliar ramp, between t = 0.06 and t = 0.0632, are 
represented. Equation (|22|) predicts index for r\ = —1.1 and index 1 for 77 = —1.0. A chaotic 
orbit of the system is also represented, and the same comment done for Figure applies here. 
Notice that for this value of E there is also a period-1 UPO, which, however, is not stabilized by 
the above values of 77. However, for 77 < —1.6, the system could choose to stabilize on either the 
period-1 or period-2 UPOs, since both are zero on the feedback control in this case. At present 
we do not have a way to predict which orbit would the system choose in a given case, or whether 
there are well defined basins of attraction. 

8 Conclusions 

We have obtained explicit, analytically computed expressions for the function g : S 1 — > C whose 
index around the origin yields the domain of control of the ETDAS method for a variety of 
orbits and feedback schemes of the buck converter. The analytical computation is possible due 
to the linear character of the differential equations of both topologies of the buck converter and 
the fact that the change of topology only introduces a delta function in the relevant equations, 
which can be translated to a delta function in the time variable. Things are not so easy for other 
DC-DC converters. For instance, for the current-controlled boost, one gets a delta function of 
the time derivative of the intensity, which is not differentiable at the topology-changing points. 
Of the three feedback schemes that we have presented, the two first are clearly superior in 
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the sense that the feedback gain needed to stabilize the UPOs increases slowly as we enter the 
chaotic zone, while the third scheme is unable to stabilize the UPOs in the chaotic regime for low 
values of r and, in any case, higher values of the feedback gain than in the other two cases are 
needed. The domain of control is quite simple in all the cases, and does not change dramatically 
as large values of r are used, although extended delay enlarges the stabilizing region in the third 
scheme. The simplicity of the domain of control for the buck converter makes it plausible that 
its general features will not change in a relevant way if nonideal elements, such as resistors in the 
inductor or in the voltage source that delivers Av are introduced, or if small noise is considered. 

We have numerically confirmed our results by simulating the time-delayed system and check- 
ing the predicted behaviour when varying the feedback gain. The simulations show that ETDAS 
is quite effective for the buck converter, although presently we don't have a way to force the 
system to choose a given orbit when several ones compatible with the imposed period are avail- 
able. 

Work is in progress to extend our results to other DC-DC converters and to include the 
effects of nonideal elements and noise. 
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